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Abstract 

The Giessen Boltzmann-Uehling-Uhlenbeck transport model with relativistic mean fields is used 
to simulate p-nucleus cohisions. Antiproton absorption cross sections and momentum distributions 
of annihilation products are calculated by varying the p coupling strength to the mean meson fields. 
Parameters of the antiproton-nucleus optical potential are extracted from the comparison of the 
model calculations with experimental data. 
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The real and imaginary parts of the antiproton optical potential are key quantities which 



determine p- nucleus scattering at low [l, l2|, |3|] and intermediate |^ energies, as well as 
the existence of deeply bound iV-nucleus states js], ^ Q, y, [o], The pN interaction in 
nuclear matter is, moreover, a challenging problem by itself, since it is subject to strong 
in-medium modifications. Indeed, in the simplest tp-approximation the real part of the p 
optical potential is repulsive and about 100 MeV high for the NN amplitude t taken at 
11| , while p-atomic phenomenology requires a strong attraction 0.0" 



threshold 



Usually the nuclear part of the p optical potential is parameterized in a Woods-Saxon 
(WS) form 



Kpt- 



Vo 



(1) 



In spite of several previous attempts to fix the p optical potential, considerable ambiguity 
still remains in its parameters. The angular distributions of elastically scattered antiprotons 
favour a shallow real part Vq = -i- 70 MeV and a deep imaginary part Wq = 70 ^ 150 



MeV in the interior of a nucleus 



mm 



15 



16l | . The Glauber and optical model 



calculations for the p absorption on nuclei [17|, Il8| also assume a negligibly small real part of 
Vopt and a strongly absoptive imaginary part. At the same time, the most recent combined 
analysis of the X-ray transitions in antiprotonic atoms and of the radiochemical data has 
produced a deep real part Vq = 110 MeV and an imaginary part Wq = 160 MeV. 

The p optical potentials from elastic scattering and p-atomic data are well determined 
only at the extreme periphery of a nucleus, corresponding to less than 10% of the central 
density 0, 9]. On the other hand, p production in proton-nucleus and nucleus-nucleus 
collisions probes the p potential deeply inside the nucleus and favours Vq = 100 ^ 200 MeV 
c„„s,ste,. w,th a dispersion relation between real and .nraginary parts of M. However, 
the microscopic transport analysis of Ref. [19( is also sensitive to rather uncertain in-medium 
elementary p-production cross sections close to threshold. 

The purpose of this work is to extract the information on a p optical potential from the 



2l| and from the data on inclusive pion 



data on p absorption cross section on nuclei 

m 

and proton production from low-energy p annihilation in nuclei |2^. Here, the absorption 
means the removal of a p from a beam caused by the annihilation, (in)elastic scattering 
and charge exchange reactions on individual nucleons. The diffractive elastic scattering on 
a nucleus as a whole is excluded from the absorption cross section. Since the absorption 



2 



requires at least one pN collision, it is sensitive to the p optical potential in a deeper region of 
a nucleus with respect to the case of diffractive elastic scattering. Indeed, the annihilation 
of 180 MeV antiprotons takes place at about half-density radius 
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24| . By the same 



argument, the proton and pion production from p annihilation on nuclei should also probe 
the p optical potential at about half-density radius. 

In our calculations, we have applied the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) 



model 



25| . This model solves a system of semiclassical kinetic equations for baryons, an- 



tibaryons and mesons coupled via collision terms and mean fields. The phase space distri- 
bution function of every particle species is projected on a set of point-like test particles. The 
coordinates Vj and space components of the kinetic four-momentum p*^ of a baryon (j = B) 
or an antibaryon (j = B) test particle are propagated in time according to the following 
Hamiltonian-like equations (c.f. 
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28l | and refs. therein): 



„*o ' 



(2) 

(3) 

with k = 1,2,3 and /i = 0, 1, 2, 3. The particles are assumed to be on the effective mass 
shell, p*^ = ^ {p^y + {rnjy. The kinetic four-momentum is defined as p*^ = pj — Vj'', where 
is a canonical four-momentum and 



^3 



g^,uj^ + gp,T'p'^ + -{B^+T^)A^ 



(4) 



is a vector field with = +1 for p and n, = — 1 for p and n, and Bb = +1, Bjj = — 1 being 
the baryon number. The field tensor in the r.h.s. of ([3]) is defined as F^^ = d^V-^ — d^V-'. 
The effective mass m* is expressed in terms of a scalar potential Sj = g^jcr as 



rrij + Sj. 



Mesonic mean fields included into the model are (/, 5') = (0,1) cu, (1,1) p , and (0,0) a 
with Qpj, and g^j being the respective coupling constants. The time component of the 
electromagnetic field, i.e. Coulomb potential, is also taken into account. The mesonic mean 
fields and Coulomb potential are calculated by solving the field equations with source terms 



provided by the currents and sca^ 



Mean Field (RMF) model 



3, 



271, 



ar density of test particles on the basis of a Relativistic 



28| . Initial positions and momenta of the test particles in 



the ground state nuclei are chosen randomly according to the spatial density distributions 
and the local Fermi momentum distribution. The proton and neutron densities are taken in 



a WS form consistent with a Skyrme Hartree-Fock systematics 



29|. 



The RMF and optical models can be related via a Schrodinger equivalent potential 



92 _ (y0)2 yO 

Re(Kpt) = S, + + -^—^ + ^i?iab, (5) 

2mj rrij 



where E'lab = yPiab + "^j ~ '"^i kinetic energy of a beam particle far away from a 

nucleus. The real part of the p optical potential becomes deeper with increasing i^iab due to 
ihe negative vector potential, in distinction to the real part of the proton optical potential 

The a-, uj- and p-nucleon coupling constants have been taken from the NL3 model 31 1 
providing a very good description of the ground states for both spherical and deformed 
nuclei. The meson-antinucleon coupling constants are quite uncertain. Following we 
introduce g^^ = -i9uN, 9pN = ^9pN, 9aN = ^9aN, where < ^ < 1 is an adjutsable 
parameter. The case of ^ = 1 corresponds to the G-parity transformed nucleon fields. For 
^ = 1, neglecting the Coulomb field, the value of the p vector potential in nuclear matter 
at the saturation density po = 0.148 fm~'^ is Vp = —308 MeV, while the scalar potential 
is Sp = —380 MeV. This gives an extremely deep real part, Re(Kpt) = —661 MeV. Below, 
we will try to find out the values of ^ which are best suited to describe p absorption and 
annihilation data on nuclei. 

The antinucleon-nucleon collision terms in kinetic equations describe the elastic scat- 
tering, inelastic production and annihilation processes: NN NN (including charge ex- 
change), NN BB + mesons, NN mesons. The cross sections of these processes are 



based on the experimental data parameterizations 32, l33|] . The NN annihilation has been 



described on the basis of a statistical model 



34 



35] . The annihilation final state includes up 



to six particles, which are various combinations of vr, 77, uo and p mesons. The annihilation 
model was originally used in Ref. for slow antiprotons, but after the proper parameter 
adjustment the model also describes successfully pion spectra and multiplicity distributions 



in pp annihilation in flight, up to pi^b — 10 GeV/c j35|]. A more detailed description of the 
NN collision channels implementation in the GiBUU model will be given elsewhere. 

To calculate collision terms, we have used a full ensemble technique of the test particle 



method (c.f. 36| and refs. therein). The full ensemble technique, in distinction to the 
cascade-like parallel ensemble one, solves the Boltzmann equation more precisely as the 
binary collisions are better localized. In the p absorption calculation on nuclei, we have 
turned off the two-body collisions of the secondary particles in order to enforce a Glauber- 




FIG. 1: (Color online) p absorption cross section on various nuclei vs the beam momentum. The 
lines marked with the value of a scaling factor show the GiBUU results. On the left three 
panels, thin solid lines represent the Glauber model calculation, Eq.Q. For the p+^^C system, 
a calculation with ^ = without annihilation is additionally shown by the dotted line on the 



lower left panel. Data are from Ref. [3| (filled boxes), Ref. [J] (filled circles), Ref. 



m (filled 



triangles) and Ref. [2l| (filled upside-down triangles). The right three panels show o"abs at lower 
jeam momenta, where the actual fit has been done. The optical model (OM) results are from Ref. 



a 



type description. 

Fig. [T] shows the p absorption cross section on ^^C, ^'^Al and ^^Cu as a function of the 
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beam momentum. The absorption cross section has been computed as 

CTabs = 27r J dbbP^^b) , (6) 


where b is an impact parameter and Pabs(&) is the probabihty of the p absorption, i.e. 
suffering at least one (in)elastic scattering or annihilation on a nucleon. The maximum 
value of the impact parameter, femax — {I-IA^^^ + 5) fm, has been chosen large enough to 
have -Pabs(&max) = in a GiBUU calculation. Without any mean field effects, Eq.([n]) can be 



reduced to a simple Glauber formula 



Glauber 
"abs 



27r J dbb[l-exp{-a,otS{b))] , (7) 



where atot is the isospin-averaged total pN cross section and 

oo 

S{b) = 2 f p{Vb^ + s^)ds (8) 







is the nuclear density integral along the straight line trajectory of a projectile with p(r) 
being the nuclear density at the radial distance r. The Glauber formula (I7l) describes quite 
well the experimental absorption cross section at very high beam momenta. However, there 
are significant deviations of Eq. (JTj) from the data at piab < 20 GeV/c. 

At piab > 0.4 GeV/c, the GiBUU calculations with ^ = are very close to the Glauber 
model, as expected. At smaller beam momenta, the attractive Coulomb potential increases 
the absorption cross section with respect to Eq. (jTj). As one can see from Fig. [H the 
agreement with the data can only be achieved when the mean meson fields are introduced. 
This can be understood as follows: In the calculation without any mean field, the beam 
particles with impact parameters larger than the nuclear radius do not experience binary 
collisions since they propagate along straight-line trajectories. Turning on the attractive 
mean field bends trajectories of the beam particles toward the nucleus. Thus, the attractive 
mean field makes a larger part of the beam flux to experience two-body collisions. 

The sensitivity of the absorption cross section to the p mean field grows with decreasing 
beam momentum. Thus we have selected the KEK data [9] at piab = 470 — 880 MeV/c to 
find the optimum value of the parameter ^ for ^^C, ^^Al and ^^Cu targets. As one can see 
from Fig. [11 ^ = 0.2 ~ 0.3 provides the best overall agreement with the data. A stronger 
attraction, i.e. larger ^, leads to an overestimation of the absorption cross section at piab < 1 



TABLE I: The scaling factor ^ of the antibaryon couphng constants and p optical potential pa- 
rameters Vq (MeV), Rniim), an (fm), Wq (MeV), Ri (fm), aj (fm) (see Eq. ([I])), obtained by 
fitting the data of Ref. [3] for different nuclei. The values per degree of freedom {F = 5) and 
the standard errors of the scaling factor ^ are given. The errors of the real depth Vb are caused by 
a variation of ^ by one standard error. The errors of all other parameters are less than 2% and are 
not shown. 



Nucleus 




XVF 


^0 




a-R 






ai 


12c 


0.22±0.03 


2.2 


153±21 


1.00 


0.63 


110 


0.97 


0.52 


27 Al 


0.21±0.04 


1.1 


162±37 


1.04 


0.64 


108 


0.99 


0.66 


64CU 


0.21±0.04 


3.3 


153±29 


1.09 


0.64 


103 


1.06 


0.65 



GeV/c. Minimizing ^2 deviation from six data points for each nucleus results in ^ values 
listed in Table [H The global fit to eighteen data points for all three nuclei produces the 
scaling factor ^ = 0.21 ± 0.03 with I F = 2.0 {F = 17). 

The quality of our calculations is visualised in Fig. [1] (right panels), where we also show 

'or a comparison the absorption cross sections from the optical model calculations of Ref. 

sl. Our absorption cross section drops with increasing beam momentum somewhat faster 
than the data do. The optical model describes the data slightly better. However, the optical 
potential of Ref. has the two free parameters, Vq and Wq, vs only one, i.e. ^ (or, 
aquivalantly, V„) in ou. n,ode,. Moreove, the fixed geometrical pa.a,nete.s of the optical 
potential of Ref. [3j] are rather arbitrary and, therefore, can be considered as free parameters 
too. 

We have also extracted the parameters of a p optical potential for the best fit values 
of ^. The real part of Vopt has been determined from Eq. ([5]) dropping the Coulomb field 
contribution in Vp. The imaginary part has been calculated as 

Im(V;pt) = -l< v,,,aZf > P , (9) 

where the averaging is done with respect to the Fermi momenta of nucleons, frei is the 
relative velocity of an incoming p and a nucleon, cr^l'^ is the total in-medium pN cross 
section computed taking into account the Pauli blocking of a final nucleon state for the 
(in)elastic scattering contribution, and p is the density of nucleons. The radial dependencies 
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of the real and imaginary parts, Eqs. andQ, have been approximated by Eq. ([T]) for 
Eiah = 0. Resulting WS parameters are listed in Table [H The mass dependence of the 
WS radii Rr, Rj appreciably deviates from a standard oc A^^^ behaviour, which is mostly 
caused by the underlying realistic neutron and proton density distributions from a Skyrme 
Hartree-Fock systematics 29|] . The extraction of the optical potential parameters is far from 
unique. In particular, the WS depths are sensitive to the assumed size parameters, and there 
is no guarantee that the WS optical potential parameters listed in Table [T] result in as good 
fit to the data as provided by the GiBUU calculation. 

One could also notice from Fig. [1] (left panels), that the BNL and Serpukhov [2^ data 
at piab = 1-6 -r- 20 GeV/c are consistent with ^ = 1, i.e. with the G-parity value of the 
real part of p optical potential. A similar result based on the data |4| has been obtained 



earlier in Ref. 



30l |. The G-parity motivated p potential is, however, not supported by more 



recent low energy data of Ref. Sf]. It would be, therefore, quite useful to perform the new 
measurements of p absorption cross sections above 1 GeV/c, which is accessible at the future 
Facility for Antiproton and Ion Research (FAIR). 

As shown in Fig. [21 further constraints on the p optical potential can be obtained from 
inclusive momentum spectra of positive pions and protons produced in p+^^C and p+^^^U 
interactions at 608 MeV/c. The two-slope structure of the pion spectra with the slope 
change at p ~ 0.3 GeV/c is caused by pion-nucleon rescattering mediated by the A(1232) 
resonance. Higher momentum pions leave the nucleus practically without interactions with 
nucleons. Lower momentum pions are either absorbed via A(1232) resonances, AA^ — > A^A^, 
or get decelerated in collisions with nucleons. 

The calculated proton spectra also change their slopes at p ~ pp, where pp = 0.27 
GeV/c is the Fermi momentum of nucleons. High momentum protons are knocked-out from 
a nucleus by energetic pions. The lower part of a proton momentum spectrum is populated 
by the slow evaporated protons produced after the fast cascading pions and nucleons have 
already left the nucleus. 

Varying parameters of the antiproton mean field influences the momentum spectra of 
annihilation products only moderately. The attraction of incoming p to a nucleus increases 
the annihilation probability. On the other hand, the invariant energy of an annihilating pN 
pair is reduced by a stronger p-attraction. As a consequence, the multiplicities and kinetic 
energies of annihilation mesons get reduced. A partial cancellation of these two effects leads 
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FIG. 2: (Color online) The angle-integrated vr"*" and proton laboratory momentum inclusive spectra 
from p interaction with ^^C and ^38^ ggg MeV/c. Calculated histograms are denoted by the 



value of the scaling factor ^. Data are from 



221. 



to a rather weak sensitivity of the momentum spectra to the p mean field. Overall, the 
calculation with ^ = 0.3 is in the best agreement with vr^ and proton momentum spectra at 
608 MeV/c. 

In conclusion, we have applied the hadron transport GiBUU model to describe p-nucleus 
interactions at the beam momenta in the range 0.4 -=-280 GeV/c. The depth of the real part 
of the p optical potential, extracted by fitting the KEK data 3| at low beam momenta is 
/n ~ 150 ± 30 MeV which is about 40% deeper than the value reported by Friedman et al 
9|. However, depths corresponding to different geometries - size parameters - may not be 
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22| favour 



directly comparable with each other. 

The annihilation spectra of positive pions and protons measured at LEAR 
even deeper real part, Vq ~ 220 ± 70 MeV. Such attractive potentials may lead to the 
cold compression effect when p penetrates deeply into the nuclear interior [7|, 



271, 



38|. 



This possibility deserves further studies in view of future FAIR experiments with antiproton 
beams. 
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FIG. 3: (addition to the PRC version of the paper) p differential elastic scattering cross section. 
The computed histograms are marked with the value of a scaling factor ^. Data are from Ref. 
Q]. With some reservations for the semiclassical nature of our model, the agreement with the 
data is fairly good within the values of the scaling factor ^ = 0.22 it 0.06, except for small-angle 
scattering at lower beam momenta, where we underpredict the experiment rather substantially. 
This deficiency is, most probably, due to neglecting the Coulomb corrections on the distant parts 
of the p trajectory (more than ~ (1.1^4^/'^ + 5) fm from the nuclear centre). The detailed study of 
elastic scattering goes, however, beyond the scope of the present work. 
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FIG. 4: (addition to the PRC version of the paper) ir'^ and proton laboratory rapidity spectra 
from antiproton interaction with ^^C and ^^^U at 608 MeV/c. Calculated histograms are denoted 
by the value of the scaling factor ^. For the 7r+ (proton) spectra, the upper lines and data points 
correspond to the transverse momenta pT > 120 (330) MeV/c, while lower lines and data points 



are computed for pT > 500 (600) MeV/c. Experimental data are from 
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